Dynamical heterogeneities close to a colloidal gel 
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Dynamical heterogeneities in a colloidal fluid close to gelation are studied by means of computer 
simulations. A clear distinction between some fast particles and the rest, slow ones, is observed, 
yielding a picture of the gel composed by two populations with different mobilities. Analyzing the 
statics and dynamics of both sets of particles, it is shown that the slow particles form a network of 
stuck particles, whereas the fast ones are able to move over long distances. Correlation functions 
show that the environment of the fast particles relaxes much faster than that of the slow ones, but at 
short times the bonds between fast particles are longer lived due to the flexibility of their structure. 
No string-like motion is observed for the fast particles, but they occupy preferential sites in the 
surface of the structure formed by the slow ones. 

PACS numbers; 82.70.Dd, 64.70.Pf, 82.70.Gg 



I. INTRODUCTION 

Colloidal gelation occurs in a system with short range 
attractive interactions, at high interaction strength for a 
wide range of density. Upon increasing the attraction, 
the system falls out of equilibrium, forming a network of 
particles with voids, that spans through the system. As 
observed by light scattering, a static low angle peak be- 
comes apparent, signaling the formation of a ramified 
structureiiS. This process has no analogue in atomic 
systems and presents a challenge in condensed matter 
physics. It is, thus, receiving increasing attention over 
the past few yeara3i4,5,6,7,8,9,io,ii,i2^i3,i4,i5^ I^s applica- 
tion in industrial processes and other sciences makes gela- 
tion even more intriguing^. Also, understanding gelation 
at low densities may lead to a full description of colloidal 
aggregation-'^^^i^ilS^iS^.. 

Colloidal gelation has been tackled recently by means 
of mode coupling theory (MCT)^^'^^, which addresses the 
glass transition both in atomic and colloidal fluida^i^iSii. 
Within MCT, gelation reduces to a regular ergodic to 
non-ergodic transition, similar to the glass transition tak- 
ing place at high density^"'-^^-^'''-^®. As predicted by MCT, 
the gel transition occurs for all densities, for higher at- 
traction strength the lower the density. Using molecular 
dynamics simulations, the predictions from MCT have 
been tested at high density, far from the percolation line, 
and good agreement has been found, showing that MCT 
correctly describes this transition2LS2422i2£. However, at 
low density this description is still under debatei^*iSi2i*S. 

Computer simulations of the glass transition at high 
density have shown also good agreement with the predic- 
tions from MCT33i34,35^ However, some differences have 
also been pointed out, especially the so-called dynamical 
heterogeneities^SiSLSi^. Particles with different mobil- 
ities are detected in the glass and in the fluid close to 
the glass, and string-like motions of particles have been 
reported. These heterogeneities in the system, which are 
more pronounced as the glass transition is approached. 



have been observed also experimentally^S*^. In this 
work, by means of computer simulations, we study if dy- 
namical heterogeneities are also present in fluid states 
close to gelation. The system under study is a colloid- 
polymer mixture, modeled by the Asakura-Oosawa in- 
teraction potential between the colloids, which has been 
shown previously to undergo a non-ergodic transition at 
high attraction strengthSiiiiiiSLffi. 

Analyzing the distribution of squared displacements, 
we show that dynamical heterogeneities indeed exist in 
the gel, and, in a simplified view, two sets of particles 
with different mobilities can be recognized. The fast 
particles have fewer neighbors than the slow ones, and 
are also smaller (a polydisperse system is simulated to 
prevent crystallization). Since the exchange of particles 
between the fast and slow populations is very slow, we 
describe the system as composed by two distinct pop- 
ulations of particles, with different mobilities. The slow 
particles form a percolating network with the same struc- 
ture as the overall system, very stable in time, whereas 
the structure formed by the fast ones relaxes faster. In 
the surface of the network formed by the slow particles, 
there exist preferential sites for the fast ones, which are 
occupied when vacant. The environment of the particles 
shows that fast particles are surrounded by other fast 
particles mainly, whereas slow ones bond to other slow 
ones. The stiffness of the structure of the slow particles 
makes it more fragile, and, at short times, bonds between 
short particles are shorter lived than those between fast 
ones. 

In this paper we analyze solely the dynamic het- 
erogeneities, and refer the reader to a previous paper 
where the MCT predictions are tested for this particu- 
lar system2L22i. In section II we present the details of 
the system and of the simulation method. Section III 
deals with the results; first we present the distribution of 
squared displacements, where the dynamic heterogeneity 
is clearly inferred. Then, the static properties of both 
populations are analyzed and their dynamics are stud- 
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ied. We then discuss the results and present the relevant 
conclusions of this work. 



II. SIMULATION DETAILS 

Molecular dynamics simulations were performed to an- 
alyze the dynamics of single particles in the system close 
to gelling. The system is comprised of 1000 soft core poly- 
disperse particles, with a short range attraction given by 
the Asakura Oosawa potential, modeling a mixture of 
colloids with non-adsorbing polymer. The core-core re- 
pulsion is given by: 



Vsc[r) 
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where ai2 — (ai -I- a2)/2, with ai and a2 the radii of 
the interacting particles. Particle sizes are distributed 
according to a flat distribution of half- width S = 0.1a, 
where a is the mean radius; ai G [0.9a, 1.1a]. The attrac- 
tion induced by the polymers, extended to take polydis- 
persity into account, reada^2iiii4: 



VAo{r) = -kBT(j)p 
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for 2(ai2 + C/5) <r< 2(ai2 
tances. Here, rji — ai/^; fj = 



^) and for larger dis- 
ivi + '/2)/2, and 0p is 
the volume fraction of the polymer. Note that the range 
of the potential is given by ^, the polymer size, and its 
strength is proportional to (f>p. For r < 2(ai2 -I- 
the potential is parabolic, connected analytically to Vao 
at r = 2(ai2 -I- ^/5), with the minimum at r = 2ai2. 
The total potential, thus, has a minimum very close to 
r = 2ai2. 

Since liquid-gas separation would prevent us from ac- 
cessing a huge area of the density-attraction strength 
plane, a long-range barrier has been added to the in- 
teraction potential. This barrier has the form: 



Vbar{r) = ksT ■ 



- 2 



r — ri 



(3) 



for < r < ri and zero otherwise. The limits of the 
barrier were set to tq — 2(ai2 + ^), and ri = 4a, and its 
height is only IfcsT'*^. The barrier raises the energy of a 
dense phase, so that liquid gas separation does not take 
place. 

The resulting total interaction potential, Vtot = Vsc + 
Vao + Vhar, which is analytical everywhere, is shown in 
Fig. n In order to show the effect of polydispersity. 
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FIG. 1: Total pair interaction potential Vtot as function of 
the radial distance r = |ri — r2| for three different particle 
pairs; a pair of particles with minimal radii ai = a2 = a — S, 
one with average radii ai = a2 — a, and one with maximal 
ai — a2 — a + S (from left to right). 



the extreme cases, as well as the average one are pre- 
sented. For small particles, not only the core interaction 
is altered, but also the attraction is weaker, whereas big 
particles interact with a stronger attraction. The max- 
imal difference between interacting pairs of particles is 
1.25 keT. 

In our simulations, lengths are measured in units of 
the average radius, a, and time in units of \fAtP^Jiv^ ^ 
where the thermal velocity v was set to \/4/3. Equa- 
tions of motion were integrated using the velocity- Verlet 
algorithm, in the canonical ensemble (constant NTV), to 
mimic the colloidal dynamics, with a time step equal to 
0.0025. Every rit time steps, the velocity of the parti- 
cles was re-scaled to assure constant temperature. No 
effect of rit was observed for well equilibrated samples. 
Equilibration of the samples was tested by the trends of 
the energy of the system and other order parametersSi, 
and by the independency of the correlation functions on 
the initial time. For the highest polymer fraction studied 
here, (/ip = 0.42, equilibration ran for 5 • 10^ time units, 
corresponding to 2 • 10^ time steps. 

The range of the attraction is set to 2^ — 0.2. The 
density of colloids is reported as volume fraction, 0c = 

Ittc"^ ^1-1- (1)^^ nc, with ric the colloid number density, 

and the attraction strength is given by the polymer vol- 
ume fraction, </>„. 



III. RESULTS AND DISCUSSION 

Dynamical heterogeneities in a system can be detected 
analyzing the distribution of the squared displacement 
of all particles at a fixed time. For homogeneous col- 
loidal fluids this distribution is single peaked, the width 
depending on the self diffusion coefficient of the particles 
and on timc"*^. However, this is not the case when popu- 
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FIG. 2: Distribution of squared displacement for <jf>p = 0.30 
(gray dotted line), <f)p = 0.40 (gray dashed line), <j)p — 0.41 
(gray solid line), — 0.415 (black dotted line), (l>p = 0.42 
(black dashed line), and (j>p = 0.425 (black solid line). In all 
cases {Sr^) — lOa^ . The thin black line is the theoretical dis- 
tribution for a system of non-interacting Brownian particles. 
Inset: Number of slow particles (open circles) and fast parti- 
cles (closed circles) in the system as a function of the polymer 
fraction. Gelation is estimated to occur at (/>„ — 0.4265. 



lations of particles with different mobilities are present. 

In Fig. 121 the distribution of squared displacements^ 
is plotted for systems at different states for increasing 
attraction strength, from time zero to time t* defined by 
{5r'^{t*)) — IQa? . At low polymer fraction, the distribu- 
tion is very similar to that of a system composed of non- 
interacting Brownian particles^. However, as the poly- 
mer concentration, 0p, is increased the distribution be- 
comes broader and two peaks are apparent, one caused by 
'slow particles' at Sr"^ ~ 5 • 10~^a^, of the order of the in- 
teraction range, 2^ — 0.2, and another one at (5r^ « 20a^ 
signaling a population of 'fast particles'. Note that the 
mean squared displacement is ((5r^) = IQo? in all cases. 
The values of t* are given in Table I for the states pre- 
sented in Fig. 121 





t* 


At 


0.30 


29.26 


2.55 


0.40 


981.6 


70.42 


0.41 


1599 


120.44 


0.415 


2162 


132.78 


0.42 


2871 


196.18 


0.425 


4911 


289.85 



Table I: Time t* , when the mean squared displacement is 



equal to lOa'^, and interval between consecutive correlators. 
At. 



The minimum in the distribution at Sr'^ — a} permits 
us to establish an unambiguous way of distinguishing be- 
tween 'fast' and 'slow' particles, instead of using an ar- 
bitrary number of fastest or slowest particles, as done in 
the study of glasses^* . The number of slow and fast par- 
ticles is plotted in the inset to the figure as a function of 
the polymer fraction. Indeed, as the attraction strength- 
ens more slow particles and fewer fast ones are present 
in the system. The gel transition was estimated to take 
place at = 0.4265 by the power-law divergence of the 
time scalers. It is interesting to note that extrapolation 
of this plot yields a noticeable fraction of fast particles 
even in the gel. 

The structure and dynamics of both types of particles, 
fast and slow, are analyzed in the following subsections. 
First, we will study the static properties, such as the dis- 
tribution of fast and slow particles, and will try to cor- 
relate the character with other properties of the particle. 
Then, we will pay attention to the motion of single par- 
ticles as well as its environment and the relaxation of the 
structures formed by the particles. Finally, we will dis- 
cuss the results, comparing them with other observations 
and theory. 



A. Static properties 

The system at 0^ = 0.40 and — 0.42, which has a 
strong distinction between fast and slow particles, is a 
heterogeneous network of particles, with voids and tun- 
nels, as presented in Fig. The distribution of the 
slow particles in the system is studied in the lower part 
of the figure, where only these particles have been shown. 
A percolating network of particles is again formed, with 
bigger voids than in the complete system. However, since 
the structure of the system is so intricate, these three di- 
mensional presentations are unclear. 

In Fig. 2] we present four slices of width 2a of the 
system, where the fast and slow particles are differently 
colored, dark grey particles are fast ones and light grey 
are slow ones. It can be seen that the slow particles form 
compact structures which are connected (in the z direc- 
tion). The fast particles, on the other hand, are also clus- 
tered but lonely fast particles are also observed, bonded 
to slow ones. The overall appearance is that no clear 
segregation between fast and slow particles is observed. 



The structure factor of the system is presented in Fig. for the same system as Figs. |21 and 0] as well as the 
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FIG. 3: Image of the system for <j!>c — 0.40 and (pp = 0.42 for 
the whole system (upper panel). Only the slow particles are 
shown in the lower panel. 



contributions from pairs of fast- fast, slow-slow and fast- 
slow particles. The partial structure factors are defined 
by: 

fc=i 1=1 

where the first sum runs over the Na fast (slow) particles, 
and the second one over the Np fast (slow) ones, and N 
is the total number of particles. The low angle peak 
observed in the total structure factor and in the fast- 
fast and slow-slow contributions is caused by the voids 
in the network of particles. In contrast, the fast-slow 
contribution, shows a minimum in this region, caused by 
the clustering of fast and slow particles. 

In the inset to this figure, the fast-fast and slow-slow 
contributions are rescaled by a factor N/Nf and N/Ng, 
the inverse fraction of fast and slow particles, respec- 
tively. Whereas the slow-slow contribution reproduces 
the neighbors peaks of the total S{q), the contribution 
from fast particles shows less pronounced maxima, indi- 
cating that there are fewer bonds between fast particles 
than average. 



Looking at the images in Fig. 01 the origin of these 
two populations can be related to the position of every 
particle. The particles in the borders of the network are 
likely to break their bonds and diffuse to longer distances, 
while those particles inside the network, are surrounded 
by many particles and cannot escape. However, it must 
also be considered that we have a polydisperse system, 
where small particles not only can go through smaller 
holes but also form weaker bonds (see Fig. 0). Both 
possibilities can be tested by correlating the displacement 
of every particle with its number of neighbors, and with 
its size. 

Fig. I^lpresents the mean number of neighbors of every 
particle averaged during the simulation up to time t* 
(upper panel) and its radius (lower panel) as a function 
of the squared displacement. Fast particles are mainly 
smaller and have fewer neighbors than the slow ones, as 
expected. However, it is also noticeable that there are 
also many small particles which are slow, indicating that 
discrimination by size has not completely occurred. On 
the other hand, many slow particles have as few bonds as 
the fast ones, showing that they are in the surface of the 
network but not breaking them and moving away. The 
mean number of neighbors of fast and slow particles are 
presented in Table II, where the distinction between fast 
and slow neighbors are made. It can be observed that 
fast particles have fewer neighbors than slow ones, and 
that the neighbors of slow particles are mainly slow ones, 
whereas fast particles are surrounded by more fast than 
slow particles. 





Total 


Fast 


Slow 


Fast 


6.07 


3.74 


2.33 


Slow 


8.37 


1.06 


7.31 



Table II: Total number of neighbors of fast and slow par- 
ticles, and number of fast and slow neighbors. 

Because there is no intrinsic difference between fast 
and slow particles, and the existence of potential fast 
particles which are slow (with a only a few neighbors), 
it is possible that an exchange of particles between these 
two populations can be taking place. For this purpose, we 
have measured 50 correlators, each one started when the 
mean squared displacement of the previous one is equal 
to a, the mean radius. In every correlator the number 
of fast and slow particles can be determined by studying 
the distribution of squared displacements at time t*, as 
shown in Fig. [21 In order to study the exchange between 
fast and slow particles, we have measured the fraction 
of fast particles which were fast in the first correlator. 
The results are presented in Fig. das a function of the 
correlator, for different polymer fractions. To make eas- 
ier interpretation of the figure the time corresponding to 
the first correlator is also given for the case (j)p = 0.42, 
and the delay time between consecutive correlators, At 
is provided in Table I. 
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FIG. 4; Slices of the system at <j!>c = 0.40 and (j) = 0.42 for z — 7.5a (upper left), z — 2.5a (upper right), z = —2.5a (lower left) 
and z — —7.5a flower right). Dark erev particles are fast ones (Sr(t*) > a^ 



and light grey means slow particles {Sr{t*) 



< a^ 




FIG. 5: Structure factor (thick grey line), and partial struc- 
ture factors for fast-fast particles (thick black line) , slow-slow 
particles (dashed black line) and fast-slow ones (thin black 
line), for (j)c = 0.40 and <j}p — 0.42. Inset: Structure factors 
rescaled to Sap{q — > oo) 1. 



Although the fraction of fast particles in every corre- 
lator that were fast in the first one decreases with time, 
the number of fast particles is constant throughout the 
simulation, as shown in the inset for (h„ = 0.42. This 



implies that indeed there is an exchange between both 
populations. However, after 50 correlators, only 30% of 
the fast particles have formed enough bonds to become 
slow, which for the (pp = 0.42 case amounts to about 100 
particles. From the point of view of the ca. 700 slow 
particles, one hundred of them have been able to break 
their bonds and difuse away. 

From one correlator to the next, only about 15 particles 
change their character, while the mean squared displace- 
ment is equal to the mean particle radius. It should be 
noted, however, that the time between correlators, Ai is 
much lower than i*, the ratio between both times being 
approximately equal to 15 in all states. Thus, Fig. 
shows that the exchange between fast and slow particles 
in time t* is below 20% of the amount of fast particles 
(about 60 particles for 0p = 0.42). 

Since this exchange between fast and slow particles 
is so slow within every correlator, both species can be 
considered as being distinguishable and stable. Thus, 
we will analyze different dynamical properties of both 
populations as a function of time. In this analysis, the 
distinction between fast and slow particles is made at 
time t* . Thus, from time i = to i*, for every correlator, 
slow particles are always slow, although at longer times 
some of them may become fast. 
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FIG. 6: Mean number of neighbors for every particle through- 
out the simulation (upper panel) and particle radius as a func- 
tion of its squared displacement (lower panel) for (j!>c = 0.40 
and <i>p = 0.42. 

B. Dynamics 

The mean squared displacement of the fast and slow 
particles is presented in Fig. |Slfor the state at </)p = 0.42, 
as well as the average over the whole system. Fast par- 
ticles are faster than slow ones at all times, even at ex- 
tremely short times, when the displacement is still within 
the interaction range. The curve for slow particles shows 
that their motion is hindered at distances about 10~^a, 
and that a subdifusive motion takes over at long times. 
In the figure we have also included the mean squared 
displacement of the 15% slowest particles, which is very 
similar to the curve of the slow particles, showing that 
the distribution of squared displacements of the slow par- 
ticles is very narrow. 

To gain more information on the motion of particles, 
we have analyzed the self part of the van Hove function, 
defined as: 

G?(r,t) = l-l^^5{r- \v,{t) - ri(0)|)^ (5) 

where, again, is the number of particles of type a in 
the system. This function provides information about the 
motion of single particles, and its deviation from Gaus- 
sian behavior is quantified by the non-Gaussian parame- 
ter, Qfa^- The self part of the van Hove function for the 



FIG. 7: Fraction of fast particles in every correlator that 
were fast in the first correlator for 0c = 0.40 and (/>p = 0.40 
(thick black line), <j!>p = 0.41 (thick black dashed line), 0p = 
0.415 (grey line), ijip = 0.42 (thin solid line) and = 0.425 
(thin dashed line) as a function of correlator (and time for 
(/f>p = 0.42). Inset: Number of particles in every correlator for 
4>a = 0.40 and (j)p = 0.42. 
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FIG. 8: Mean squared displacement of the system (solid black 
line), considering only the fast particles (dashed black line), 
slow particles (solid grey line) and the 15% slowest particles 
(dashed grey line). 



system at 0c = 0.40 and (f>p — 0.42 is presented in Fig. |5| 
(upper panel) for exponentially increasing times. 

The van Hove function for the system shows the typi- 
cal behavior of close-to-glass systems; at short and inter- 
mediate times, the van Hove function has a maximum at 
the estimated localization length, of the order of the cage 
size, or the interaction range in our case. At long times, 
the particles break free of their cages (network of bonds) 
and diffuse to longer distances. However, when plotted 
in logarithmic scale, see the inset to upper panel, it can 
be observed that a non-negligible tail appears at long 
distances. This feature is caused by particles moving at 
long distances even at short times, i.e. the fast particles. 
Thus, the van Hove function has been calculated consid- 
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FIG. 9: Self part of the van Hove function for 0c — 0.40 and 
(/f>p = 0.42 (upper panel) for exponentially increasing times, 
t„ = 0.25 • 2^^"~^\ for the fast particles only (intermediate 
panel) and slow particles only (lower panel). Insets: The same 
as in the main graph in logarithmic scale. 



ering fast particles only (intermediate panel), and slow 
particles only (lower panel). It is worth remembering 
that the distinction between fast and slow particles was 
made when the mean squared displacement was equal to 
lOa^, which for this state corresponds to t* — 2871. 

Indeed, the fast particles are responsible for the tail at 
long distances, but also fewer fast particles are trapped 
in the network of bonds (see the height of the peak) 
and those which are imprisoned escape earlier than av- 
erage. On the other hand, there are more slow parti- 
cles confined, and for much longer times, and their van 
Hove function does not present the long distance tail. 
It is interesting to observe that at long times, there are 
just a few fast particles still close to their original posi- 
tions and the van Hove function is practically flat (note 
the dashed line in the second panel close to the x-axis), 
whereas most of the slow particles are still at a distance 
of the order of the interaction range. The latest time 
presented, t = 4096, is after the distinction between fast 
and slow particles was made, and therefore, some ex- 
change between the fast and slow populations may have 
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FIG. 10: Bond correlation function for (j)c — 0.40 and 
(l)p = 0.42 for the whole system (solid black line), for 
bonds between fast particles (dashed black line), bonds be- 
tween slow particles (solid grey line) and bonds between fast 
and slow particles (dashed grey line). Inset: Bond corre- 
lation function for the systems at (fic — 0.40 and (j)p — 
0.30, 0.40, 0.41, 0.415, 0.42, 0.425 (from left to right). 



taken place. This is the origin of the tail at long dis- 
tances observed in the function for the slow particles at 
this time. 

Since the driving process for the gel transition is bond 
formation, we study now the differences between fast and 
slow particles from the point of view of their bonds and of 
their environment. The bond correlation function, ^sit), 
is defined as the fraction of bonds that have uninterrupt- 
edly existed from time t = 0. This function is presented 
in Fig. ^|for different states approaching the gel tran- 
sition (inset to the figure). As the polymer fraction in- 
creases, the bonds between particles are more stable, not 
only because of the increasing attraction strength, linear 
with (/)p, but also due to the higher number of neighbors, 
forming an intricate network. This collective character 
becomes apparent also since this function departs from 
the simple exponential more and more. 

In the main body of Fig. ^|the bond correlation func- 
tion for (j)p = 0.42 is analyzed in more detail, presenting 
the contributions from bonds between two fast particles, 
two slow ones and bonds between fast and slow particles. 
The fast-fast bonds break faster at intermediate and long 
times, showing that these particles move as single parti- 
cles, and not in clusters of fast particles. In contrast, the 
most stable bonds are those formed between two slow 
particles, since they hardly move. However, this trend is 
inverted at short times; there, the bonds between slow 
particles are the easiest to break, while those between 
fast ones are more stable. This feature can be rational- 
ized considering that slow particles live in a rigid envi- 
ronment, where the motion of a single particle, even if 
small, can break bonds, as the others are not able to 
reaccommodate to keep their bonds alive. For fast parti- 
cles, the loose environment allows some fluctuation, and 
their bonds are more stable at these short times. 
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FIG. 11: Environment correlation function for (j)c = 0.40 and 
<pp = 0.42. The lines have the same meaning as Fig. 1101 Inset: 
Environment correlation function for the states as Fig. IIUI 

Although the bonds between slow particles are the 
most fragile at short times, the environment of these par- 
ticles does not change, i.e. their likely neighbors are the 
same for very long periods of time. This fact can be 
recognized by studying the fraction of neighbors of the 
particle that were neighbors a,t t — 0. This environment 
correlation function, $e(i) is plotted in Fig. ^2 for dif- 
ferent states (inset to the figure) and in more detail for 
the (j)p = 0.42 case. Because this function measures the 
relaxation of the environment of every particle, instead 
of the lifetime of the bonds, it decays much slower than 
$B, as observed in comparing Figs. 1101 and [TTI 

As expected, the environment correlation function de- 
cays slower as the attraction strengthens, and in the slow- 
est case, (j)p ~ 0.425, it has decayed to about $e = 



0.5 only, although the mean squared displacement was 
((5r^) = 20a^. The contribution to the global behavior 
from fast neighbors of fast particles decays much faster 
than average, as that of slow- fast pairs (or vice versa). In 
contrast, the neighborhood of the slow particles formed 
by other slow particles is hardly altered, keeping more 
than 95% of their neighbours in average. It is impor- 
tant to note that although the bonds between slow parti- 
cles can break, their neighbours do not change. The sys- 
tem, thus, is constantly rearranging, breaking and form- 
ing new bonds, but the overall structure is quite stable. 

The behavior of the correlation functions <i>s(i) and 
$e(i) for the bonds and environment of the fast parti- 
cles, indicates that these particles break their bonds with 
their neighbors and change them, showing that their fast 
motion does not take place inside fast clusters. Neverthe- 
less, since the process of breaking single bonds is much 
faster than neighborhood relaxation, the restructuration 
of bonds alluded above also takes place for these parti- 
cles. In order to understand in more detail the behavior 
of these particles, we have analyzed the distinct part of 
the van Hove correlation function, defined as: 



^^"^"''^ = NjTp (^^^5(r- IrKi) -rj(O)l)^ (6) 

where N, Na and Np have the same meaning as above. 
For a — f3, the fraction in front of the summation reads 
N/Na{Na — 1), and the case « = j is omitted in the 
summation. For long times, this function tends to one, 
G{t — > oo) — > 1, as the structure of the system at t = 
relaxes. 



In Fig. I12l we present this function for the same times 
as Fig. El for the whole system as well as the contribu- 
tions from the different combinations. As expected from 
the discussion above, the function for the whole system 
does not reach its limit at long times, since the struc- 
ture does not relax in the time of the simulation. It is 
also interesting to note that the function shows a peak 
at the origin, whose height grows with time. This peak 
is caused by particles occupying exactly the same place 
where other particles were at previous times. In glasses, 
this feature is recognized as a signature of hopping, or 
string-like motions^. 

The origin of the peak at r = can be understood 
studying the contributions from the different pairs of par- 
ticles. The peak appears only in the fast-fast contribu- 
tions, though it is present at times longer than t* in the 
other contributions. It can be rationalized as fast parti- 
cles occupying the place where other fast particles laid. 



but string-like motions are not observed (by direct vi- 
sual observation of the simulation). The mean number 
of neighbors of these particles while occupying the former 
place of another fast particles is 6.54, with 3.07 fast neigh- 
bors and 3.48 slow ones. The ratio of slow-fast neighbors 
of these particles is much larger than the average over all 
of the fast particles (see table I), showing that they are 
much closer to (in contact with) the structure formed by 
the slow particles. We thus see preferential sites in the 
boundary of this network, where fast particles can attach. 
However, at long times, the function shows the expected 
limit, G{t — » Qo) 1, except for the peak at the origin, 
showing that the structure formed by the fast particles 
has disappeared. 

The distinct part of the van Hove function from slow- 
slow pairs of particles shows little evolution with time, as 
expected, since the structure formed by these particles 
is quite rigid and stable. For t = 4096, t > t* , a, low 
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peak at the origin appears, indicating that hopping can 
be detected even for the sfow particles for fong times. 
The slow particles that change to fast ones, or vice versa, 
cause the peak at r = in the contribution from the fast- 
slow pairs. The neighbor peak is in this case lower than 
in the other contributions, since there are few contacts 
between fast and slow particles. 

Finally, we will analyze the relaxation of the structures 
formed by the fast and slow particles. Figs. [TUI [TTl and 
^|have already provided information about this feature, 
though studying it locally. The relaxation of a struc- 
ture is better observed by the (coherent) density-density 
correlation function: 

<^-P{q,t) = _ /^^exp{zq. [ri(i) - rj(0)]}) (7) 

where the brackets indicate ensemble averaging, and q 
is the wave vector, giving the typical size of the struc- 
ture. The self or incoherent part of this function, $"(t) 
is calculated taking i = j. Both the coherent and inco- 
herent correlation functions are plotted in Fig. ^|for the 
system at 4>c = 0.40 and 4>p = 0.42 at the wavevectors 



of the first peak in the structure factor, qa = 1, (upper 
panel) and at the neighbor peak, qa = 3.9 (lower panel). 
The former value has been taken as representative of the 
macro-structure of the colloidal gel, whereas the latter 
gives again information about the local environment of 
the particles. 

The structure formed by the slow particles indeed 
changes very little in the time of the simulation, as ob- 
served in their contribution to the coherent part of the 
intermediate scattering function at qa = 1, which de- 
cays less than a 5% of its initial value. In contrast, the 
contribution from fast particles, decays much faster, as 
expected from the results above. It is interesting to com- 
pare the coherent and the incoherent parts of the func- 
tion. At this wave vector, the self parts decay a decade 
faster than the coherent parts. For the higher wave vec- 
tor, qa — 3.9, the correlation functions decay faster than 
at the previous q, implying that the structures at the cor- 
responding length relax faster than the overall colloidal 
structure. Again, the contributions from slow and fast 
particles show the expected behaviour, although the re- 
laxations are faster in this case. 
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FIG. 13: Coherent (thick lines) and incoherent (thin lines) 
density correlation function for (j)c = 0.40 and 0p = 0.42, at 
qa = 1 (upper panel) and qa = 3.9 (lower panel). Contribu- 
tion from the slow particles (dashed lines) and from the fast 
ones (dotted lines) are presented. 



C. Discussion 

The results presented in the previous subsections lead 
to a dynamically heterogeneous picture of a colloidal gel. 
In a simple description, it can be considered as composed 
by two distinct populations of particles: 'fast' and 'slow'. 
Although there is no intrinsic difference between both 
particles, the interchange between them has a very low 
rate. The slow particles form a rigid network of particles, 
very stable with time; the others, quite movable, move 
using sites where other fast particles laid, since the struc- 
ture formed by these particles relaxes ten times slower 
than the mean time for particle relaxation. As the attrac- 
tion strengthens, i.e. the gel transition is approached, the 
distinction between fast and slow particles becomes more 
pronounced, and the fraction of fast particles decreases, 
as observed at fixed mean squared displacement. 

These results were tested by analyzing the evolution of 
the squared displacement distribution up to longer times 
(10'* time units). Only a few correlators were consid- 
ered. For the system at (fip — 0.42, the peak of fast 
particles moves continuously to longer distances reach- 
ing the maximum distance in the simulation (one half of 
the box edge ^ 120a^). The peak of the slow particles 
slowly moves to higher distances, but stays of the order of 



10~"'^a^. This distinction up to long times indicates that 
the structure of the slow particles does not relax in the 
simulation time-window, as observed in the intermediate 
scattering function, Fig. [T^ 

This fact may raise the point of how well equilibrated 
our systems are. The slowest system, 0p = 0.42, was 
equilibrated for 5 • 10^ time units, which is more than ten 
times t* for this state. The correlation functions mea- 
sured at different times, the ratio of fast/slow particles, 
as well as other quantities, did not show any dependency 
on the initial time (between 5 ■ 10'* and 6 ■ 10^). Ap- 
parently, this state is, thus, well equilibrated. For lower 
polymer fractions, (pp < 0.41, the whole system relaxes 
within the simulation time-window. We cannot ascer- 
tain, but expect that the coherent and incoherent corre- 
lation functions of Fig. ^| decay via the exchange of fast 
and slow particles shown in Fig. [7| If so, these corre- 
lation functions should decay eventually to zero, so that 
the systems we have studied remain on the liquid rather 
than the gel side of the nonergodicity transition, as as- 
sumed in the analysis we have made in this paper, and 
previously"^''. However, since the contributions from slow 
particles to the correlation functions decay only a little 
within the observable time window, we cannot rule out 
a scenario where the correlation functions decay not to 
zero, but to a finite value, at timescales beyond those 
studied here. Such a plateau could arise in the presence 
of a persistent structural component comprised of a sub- 
set of the slow particles. There is no way to rule out such 
a scenario other than by extending the simulation runs 
by at least one further decade in time, which is beyond 
the scope of this paper. 

It is interesting to note, nevertheless, that the overall 
behavior of the system is correctly described by MCT— . 
The localization length obtained in the simulations and 
in the theory are in agreement, showing that bond forma- 
tion is the leading process in colloidal gelation. The be- 
havior of the self part of the density correlation functions 
obey the results obtained in MCT, both in time-scaling, 
and the specific form of the decay. The behaviors of the 
time scale of the correlation function as a function of 
the attraction strength and of the wavevector, also fit 
in the predictions from MCT. However, MCT gives in- 
formation through correlation functionsSiiSSi^, and thus 
cannot give insights into the origins of the dynamical het- 
erogeneities presented here. 

The fast particles fluidize the system, and allow for re- 
structuring of parts of the system, causing a decay in the 
correlation function, which can be described using MCT. 
As the gel transition is approached, there are more and 
more slow particles, forming the rigid network, but from 
our observations, we cannot say what are the limits of the 
fractions of fast or slow particles at gelation. Whereas 
MCT is able to describe the behavior of the system close 
to the gel transition, though only for averaged quanti- 
ties, it is expected that the dynamical heterogeneities 
described here will become important closer to the tran- 
sition, and deviations from the MCT behavior should be 
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expected. 

The strong dynamical heterogeneities presented here 
appear to be caused by the structural heterogeneities in 
the system. Because the latter are absent in the glass 
transition driven by steric hindrance, we do not expect 
such big effects in this transition. The string-like mo- 
tions reported in, e.g. Lenard Jones systems, which also 
fluidize the system, are caused also by structural het- 
erogeneities, although of the size of a single particle. 
Thus, their effects in the behavior of the system is not 
as dramatic as in the gel case. Presumably, the effects of 
the two populations will be more important for gelation 
at lower packing fractions, where the structural hetero- 
geneities are also more important. 

IV. CONCLUSIONS 

We have studied in this paper the dynamics of a col- 
loidal gel, showing that its microscopic behavior is far 
from homogeneous. Some particles are faster than oth- 
ers, and two populations can be discerned, according to 
their squared displacement at long times. Fast particles 
have fewer neighbors than slow ones, and thus, can es- 
cape their bonds easier. Because there is no structural 
difference between them, some fast particles become slow, 
when trapped by others, and some slow ones finally break 
free, and become fast. However, the exchange rate be- 
tween fast and slow particles is very low and allows an- 
alyzing the system as composed of two different popula- 



tions. Indeed, this oversimplified two population picture 
gives an adequate description of the main features of the 
system. 

The motion of the fast particles does not occur in clus- 
ters or arms rearranging, but they move singly, as ob- 
served in the bond correlation function. Their bonds 
are the weakest at intermediate and long times, and the 
structures formed by the fast particles relax faster than 
average. In contrast, in the network formed by the slow 
particles the bonds between slow particles are very frag- 
ile at short particles due to the stiffness of the network. 
The intermediate scattering function at the wave vector 
related to the structure of the gel, shows that this struc- 
ture is very stable with time, while the structure formed 
by the fast particles relaxes much quicker. 

These observations on the microscopic dynamics of the 
gel cannot be accounted for in MCT in its present form, 
although the average behavior is correctly described by 
this theory. Close to the gel transition, we expect, the ef- 
fects caused by the distinct populations would dominate 
the behavior of the system and the MCT picture would 
break down. 
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